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PREFACE 

This  Memorandum  is  the  result  of  further  study  of  the  trajectories 
of  glide  reentry  vehicles.  Related  Memoranda,  wiiich  were  based  on 
previous  investigations  of  the  use  of  the  atmospheric  forces  to  effect 
changes  in  the  orbital  plane  for  winged  spacecraft,  include  RM-4391-PR, 
Minimum  Energy  Loss  Heading  Changes  for  Hypersonic  Flight  from  Orbit, 
by  Russell  D.  Shaver  and  RM-3231-PR,  The  Synergetic  Plane  Change  for 
Orbiting  Spacecraft .  by  F.  S.  Nyland. 

This  present  study  ic  part  of  several  continuing  analyses  related 
to  space  operations.  It  should  be  of  interest  to  planners,  analysts  in 
the  fields  of  atmospheric  reentry,  designers  of  advanced  space  vehicles, 
and  persons  concerned  with  maneuvering  space  vehicles  and  their  utility. 
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SUMMARY 

This  Memorandum  presents  the  elements  of  a  model  for  predicting 
the  turning  trajectories  of  hypersonic  reentry  gliders  by  approximate 
closed  form  solutions.  Vehicle  control  Is  assumed  to  be  exercised  by 
varying  the  angle  of  attack  and  bank  angle.  In  this  Investigation  the 
control  variables  are  assumed  to  be  constant  for  each  particular  tra¬ 
jectory,  but  the  lift-to-drag  ratio  (L/D)  has  been  varied  from  1.0  to 
4.0  and  the  bank  angle  has  been  permitted  to  become  as  large  as  90  deg. 
Approximate  closed  form  solutions  are  compared  with  earlier  results  of 
a  study  of  minimum  energy  hypersonic  turns. 


vli- 


CONTENTS 

PREFACE .  Ill 

SUMMARY .  V 

SYWOLS  .  lx 

Section 

I.  INTRODUCTION .  1 

II.  EQUATIONS  OF  MOTION .  3 

III.  TRAJECTORIES  FOR  BANK  ANGLES  LESS  THAN  90  DEGREES .  6 

IV.  TRAJECTORIES  FOR  BANK  ANGLES  OF  90  DEGREES .  22 

V.  HEATING  TRENDS  .  35 

VI.  FINAL  coKMErrrs  .  37 


REFERENCES 


43 


ix- 


a 


A 

b 


D 

g 

h 

L 

n 

R 

o 

s 

t 

u 

u  , 
1 

u 

o 

u 

r 

W 

X 

y 

z 


SYMBOLS 


2 _ 

(L/D)  sin  4 

aerodynamic  reference  area 

8  Pq  CpA 
2  B  W  sin  (“y) 

drag  coefficient 


lift  coefficient 


drag  force 

acceleration  due  to  gravity 

altitude 

lift  force 


integer 

earth  radius  (3442  n  mi) 


path  length 
time 

vehicle  velocity 
final  velocity 

initial  velocity 

satellite  velocity  at  200,000  ft  (25,815  ft/sec) 

reentry  velocity 

vehic le  weight 

downrange  distance 

siderange  distance 

logg  *  a  dummy  variable 
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6  ■  inverse  of  atmospheric  scale  height  (1/24,000  ft) 

V  >  flight  path  angle  (measured  positive  upwards) 

a  ■  e  ■  density  factor 

■  sea  level  density  of  air  (0.07648/32.2  slugs/ft  ) 

0  ■  bank  angle 

u>  ■■  heading  angle 


I.  INTRODUCTION 


In  past  stuoies,  several  authors  have  presented  closed  form 

solutions  to  describe  the  turning  trajectories  of  lifting  vehicles 

(1  2  3) 

entering  the  atmosphere.  *  *  Loh  achieved  approximate  closed 

form  solutions  by  specifying  the  shape  of  the  ground  track,  a  minor 
circle.  Slye  showed  approximate  solutions  for  constant  bank  angle 
control  of  vehicles  with  low  lift-to-drag  ratios  (0  <  L/D  <  2) ,  but 
his  main  emphasis  was  on  maximizing  the  side  range,  rather  than  mini* 
mixing  energy  loss  for  a  given  heading  change.  In  this  Memorandum, 
approximate  closed  form  solutions  will  be  developed  that  may  be  applied  to 
a  greater  range  of  lift*to*drag  ratios  where  the  vehicle  bank  angles 
of  interest  will  be  particularly  steep,  even  up  to  90  deg.  The 
approach  here  is  quite  similar  to  the  one  used  by  Slye,  but  it  additionally 
considers  near-ballistic  reentries  in  which  all  lift  is  used  to  change 
the  v'^hlcle  heading. 

The  Interest  in  constant  bank  angle  control  is  basically  two* 

fold.  Such  a  control  law  is  easy  to  mechanize  by  either  automatic  or 

pilot  control.  Secondly,  a  recent  investigation  of  turning  trajectories 
(4) 

by  Shaver  showed  that  such  a  control  law  may  be  nearly  optimal  for 
synergetic  plane  changing^^^  where  the  optimal  path  is  defined  as  that 
trajectory  which  minimizes  the  energy  loss  during  the  aerodynamic 
portion  of  flight  for  a  given  vehicle  heading  change.  Shaver's  work 
also  indicates  that  the  optimal  control  law  for  the  angle  of  attack 
is  for  the  vehicle  to  fly  at  an  attitude  that  results  in  the  maximum 
lift*to*drag  ratio.  Thus,  as  a  result  of  this  recent  work,  the 
problem  is  no  lunger  determining  optimal  control  laws,  but  rather 
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finding  a  simple  way  of  describing  flight  trajectories  that  minimize 
energy  losses*  This  Memorandum  derives  the  descriptive  material 
necessary  for  making  preliminary  determinations  of  the  trajectory 
shapes.  A  constant  bank  angle  control  scheme  that  assumes  a  constant 
lift-to-drag  ratio  is  used. 

In  addition  to  synergetic  plane  changing,  there  are  other  appli¬ 
cations.  There  has  always  been  an  interest  in  the  side  range  that  a 
reentering  vehicle  can  achieve,  and  such  calculations  will  be  presented 
later  in  this  Memorandum. 
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II.  EQUATIONS  OF  MOTION 


The  equations  of  motion  used  in  this  analysis  have  been  written 
as  a  resolution  of  accelerating  influences  along  the  flight  path, 
along  a  line  perpendicular  to  the  flight  path  in  a  generally  up  and 
down  direction,  and  in  a  direction  out  of  the  flight  plane.  The 
accelerations  in  these  various  directions  are  assumed  to  be  control¬ 
lable  by  a  variation  in  the  bank  angle  and  the  lift* to-drag  ratio. 

These  equations  also  include  the  assumption  that  the  acceleration  due 

to  gravity  does  not  change  very  nuch  in  the  altitude  regime  of  interest. 

2 

Therefore,  the  orbital  velocity  may  be  written  as  u^  •  where 

is  the  sum  of  the  earth*8  radius  and  reference  altitude,  and  g  Is 
the  acceleration  due  to  gravity  at  the  reference  altitude. 

The  accelerations  along  the  flight  path  may  be  expressed  as 


g  dt 


-  sin 


Y 


w 


(1) 


where  u  is  the  vehicle  speed,  t  is  time,  y  is  the  flight  path  angle 
in  the  vertical  direction,  D  is  the  drag  force,  and  W  is  the  vehicle 
weight.  Fig.  1  shows  the  flight  geometry. 

Along  the  local  vertical, 

g  ^  “  W  ^ 

where  L  is  the  lift  force,  4  is  the  bank  angle,  and  u^  is  the  orbit 
speed  at  a  reference  altitude  of  200,000  ft. 

Crossways  to  the  flight  path,  the  heading  change  with  time  is 


-Earth's  limb 


Fig.  1— Flight  trajectory  coordinates 
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u  ^ 
g  dt 


0 


(3) 


where  cu  Is  the  flight  path  angle  sideways  (heading  angle)  of  the 
vehicle.  Equation  3  is  a  balance  of  forces  in  the  plane  established 
by  the  center  of  curvature  and  two  points  very  close  to  each  other  on 
the  trajectory,  but  does  not  include  the  higher  order  terms  introduced 
by  the  earth's  curvature. 
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III.  TRAJECTORIgS  FOR  BANK  ANGIZS  IZSS  THAN  90  DEGREES 

In  this  section,  solutions  to  the  equations  of  motion  will  be 
found  for  turning  trajectories  where  the  bank  angle  is  less  than  90 
deg.  Whenever  possible,  closed  form  solutions  will  be  sought  if  only 
minor  errors  are  committed  in  making  approximations.  For  one  part  of 
the  following  trajectory  description,  such  a  solution  is  apparently 
not  possible;  the  results  of  numerical  integration  will  be  presented 
for  a  large  range  of  parameter  variations. 

Our  assumptions  are  quite  similar  to  those  commonly  used  in 
analyses  of  lifting  trajectories.  It  is  assumed  that: 

1)  the  flight  path  angle  (hereafter  referring  only  to  the 
up  and  down  direction)  is  very  '^mall 

2)  the  acceleration  due  to  gravity  changes  only  slightly  in 
the  regions  of  interest,  and  that  such  changes  may  be  neg¬ 
lected 

3)  the  predominant  deceleration  influence  is  due  to  drag 

4)  the  flight  path  angle  does  not  vary  greatly  (d  y/dt  “  0) 

5)  the  control  variables,  L/D  and  (J,  are  constant  throughout 
the  flight  regime. 

These  assumptions  appear  to  introduce  only  small  errors  in  the 
solutions  to  the  equations  of  motion  provided  there  is  some  component 
of  lift  along  the  local  vertical,  so  that  the  path  angle  can  be  reasonably 
maintained  at  small  values  (y  <<  D/W) .  For  the  particular  case  in  which 
there  is  no  component  of  lift  along  the  local  vertical,  one  must  make 
other  assumptions  to  find  closed  form  solutions.  Thus,  a  later  part  of 
this  Memorandum  will  be  devoted  to  solving  the  equations  of  motion  when 
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the  bank  angle  is  90  deg  and  there  is  no  component  of  lift  along  the 
local  vertical. 

Under  these  assumptions,  the  equations  of  motion  become 

i  ^  .  2.  .  _1 _ 

g  dt  W  2  g  ds 

L  A  ^  2,2 

77  cos  0  »  1  -  u  /u 
w  o 


(la) 

(2a) 


g  dt 


L 

W 


sin  d 


(3«) 


Combining  equations  la  and  3a  to  eliminate  time  as  a  variable 
and  integrating  from  an  initial  velocity  where  the  turn  is  started, 

0)  -  ^  sin  ^  log^  (uj^/u)  (4) 

This  result  indicates  that  the  turning  rate  is  highest  at  the  end  of 

a  turn  or  that  the  turning  rate  increases  as  the  vehicle  slows  down 

(Fig.  2).  The  largest  heading  changes  can  be  accomplished  with  maximum 

L/D  and  very  steep  bank  angles.  From  equation  2a,  the  normal  load 

factor  also  Increases  as  the  vehicle  slows  down.  These  trends  are  the 

(4) 

same  as  those  noted  by  Shaver  In  optimizing  synergetic  plane  changes. 

To  determine  »'he  altitude-velocity  profile,  the  lift  force  may 
be  rewritten  as 


^  “  2  Po 


9^  2  .  . 

u  Cj^A 


3 

where  is  the  sea  level  air  density  (slugs/ft  )  ,  0  is  the  inverse 
of  scale  height  of  the  atmosphere,  h  is  the  altitude,  is  the  lift 
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Fig. 2— Heading  angle-velocity  relation 
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coefficlent,  and  A  it  the  vehicle  reference  area.  Substituting  this 
statement  into  equation  2a,  then 


3h 

e 


2/  2 

(9  (  .  9,0} 


1 


u^/u^ 

o 


— )  COS  0 


Solving  for  3h,  and  separating  Into  convenient  quantities,  note  that 


0h 


W 

i  Po  ■  e  Cpt 

o  *- 


-  log^  (1/C08  0)  (5) 

From  this  result  It  Is  possible  to  construct  a  simple  calculating  tool 
that  can  be  used  to  find  the  alti tude -velocity  profile  for  constant 
bank  angle  control  routines  (Fig.  3).  This  figure  shows  contours  for 
various  bank  angles  with  the  lift  parameter,  W/C.A,  set  to  unity.  For 

la 

other  values  of  W/C.A,  a  scale  Is  given  that  shows  the  altitude  decre- 

la 

ment  for  variations  In  this  parameter.  Also  shown  are  flat  plate 
temperatures  to  Indicate  the  heating  trend,  should  this  be  of  Interest. 
Note  here  that  because  the  fundamental  shape  of  the  equilibrium  glide 
profile  Is  unchanged  by  variations  in  the  bank  angle,  the  nAxlmum 
heating  rate  for  constant  bank  angle  trajectories  appears  to  be  less 
than  for  minor  circle  turn  trajectories.  Minor  circle  profiles  are 
such  that  the  vehicle  will  be  at  lover  altitudes  than  would  be  the  case 
for  a  constant  bank  angle  contour  at  the  same  velocity  In  the  high 
velocity  regime  where  the  maximum  heating  rate  occurs.  This  assumes 
that  the  W/C^A  the  same  in  each  case.  Also,  the  velocity  of  max¬ 
imum  heating  rate  is  shifted  to  higher  values  for  the  minor  circle 
tum.^^^  Such  Is  not  the  case  for  constant  bank  angle  control,  and 
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450. 


n;,^  I  I  I  m  I  m  I  I  I  I  1  I  i  t 

Th€  curves  shown  here  apply  to  gliders  with  W/CiA=  I  Ib/ft. 
'For  other  values  of  this  porometer,  grophicolly  subtroct  the 
-altitude  correction  indicoted  by  the  ouxiliory  scale.  A  similar 
device  for  minor  circle  turns  will  be  found  in  Ref.  5 


0.4  0.6 

Normalized  velocity  (u/uq) 

Fig. 3— Glide  reentry  trajectories  for  constant  bank  angle  control 
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the  point  of  roaxlnum  heating  would  still  be  the  same  as  determined  in 
Gaz ley's  analysis. 

Vehicle  time  of  flight  can  be  found  by  writing  an  expression  for 
the  drag  force, 

„  1  -0h  2  ^  ^ 

D  ■  T  p  e  u  C_A 
2  o  D 

fih 

and,  substituting  the  expression  for  e  from  above,  the  differential 
equation  relating  velocity  and  time  of  flight  is  generated. 


A. 

g  dt 


1  -  u^/u 


L  A 
^  cos  0 


2 

o 


Integrating  from  an  Initial  velocity  and  defining  the  time  at  the 
initial  velocity  as  t  ■■  0,  the  time  of  flight  is  given  by 

u  (L/D)  cos  0  (1  •  u/u  )  (1  +  u  /u  ) 

°  o  i  o 

Selected  times  of  flight  for  various  initial  velocities  are  shown  in 
Fig.  4.  Note  from  this  result  that  the  time  decreases  to  zero  when 
the  bank  angle  approaches  90  deg  or  L/O  approaches  zero.  Under  actual 
flight  conditions,  this  trend  would  not  hold.  The  equations  yield  such 
a  result  because  the  path  angle  would  not  be  small  when  (L/D)  cos  0 
approaches  zero;  this  violates  an  assumption  msde  earlier.  Under  these 
conditions,  the  trajectory  would  more  closely  reseinble  a  ballistic 
entry  into  the  atmosphere  than  a  lifting  trajectory.  Thus,  this 
analysis  would  probably  be  in  error  when  vehicles  with  low  L/D  are 
very  steeply  banked;  that  is,  when  the  component  of  lift  along  the 
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local  vertical  is  small.  A  separate  analysis  will  be  performed  in  a  later 
section  of  this  paper  to  consider  such  flight  regimes.  In  this  section, 
the  bank  angle  will  be  varied,  but  only  up  to  a  value  slightly  less  than 
90  deg  for  the  range  of  lift-to-drag  ratios  of  interest  (1  ^  L/D  ^  4). 

The  solution  for  the  path  length  is  very  similar  to  that  for  the 
time  of  flight,  in  that 

ds  (L/d)  ^08  0 

2 

Integrating,  and  noting  u^  ■  R^gi 


1  2,  2 

a  ■  2  R  (L/D)  cos  0  log  -  7~  2 -  O) 

®  1  -  u .  /u 

1  o 


The  path  length,  s,  is  shown  for  various  parameter  values  in  Fig.  3* 

Note  that  for  initial  velocities  of  about  0.95  u  ,  most  of  the  distance 

o 

to  be  flown  has  been  covered  by  the  time  the  vehicle  slows  down  to 
u  B  0.6  u^.  Thus,  it  would  appear  that  the  effects  of  L/D  for  covering 
distance  is  greatest  at  high  velocities,  and  thus  a  high  L/D  at  low 
velocities  does  not  contribute  nuch  in  large  scale  maneuvers. 

The  downrange  and  siderange  distances  for  constant  bank  angle 
control  are  found  by  noting  that 


dx 

du 


cos  U), 
du 


du 


du 


sin  ub, 


where  x  is  the  distance  traveled  in  the  original  direction  of  morion 
(downrange  distance)  and  y  is  the  distance  traveled  perpendicular  to 


L/D  cos 
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the  original  direction  of  motion  (sidcrangc  distance).  Differentiating 

y 

equation  7  with  respect  to  velocity. 


g 


du 


u 


g/p) 


cos 


1  - 


2/  2 
u  /u 

o 


and  substituting  into  the  differential  equations  for  crossrange  and 
downrange  distance  (to  eliminate  the  variable  of  vehicle  heading  angle), 
a  pair  of  differential  equations  is  generated  relating  velocity  and 
distance. 

g  ^  — 2~  ^  ^  ^°®e  “i^“^ 

1  -  u  /u 

o 

g  ^  *  - ^^2 — 2  ^  ^  ^°^e  “i^“^ 

1  -  u  /u 

o 

These  equations  may  be  solved  by  making  the  substitutions 

z  •  k  log  (u/uj^),  k  «  “  sin 

and  expanding  the  result  into  an  infinite  series  that  is  uniformly  con¬ 
vergent  in  the  region  of  interest.  After  integrating  term  by  term  and 
imposing  limits, 

00 

g  X  “  ^  ^(u/u  )^"  sin  sin  0  log  —  +  tan  ^  a  n) 

K  O  U  U  , 

O  X 


E 

n*l 


2n 


a  n 


^  1  +  a^  n^ 


(10) 
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2 

where  a  "■  - — .  " 

(L/D)  sin  ^ 

These  functions  are  shown  in  Figs.  6,  7,  8,  and  9  for  L/Ds  ranging  from 

1  to  4.  Velocity  is  shown  as  a  parameter,  with  the  initial  velocity 

assumed  to  be  0.98  u^.  The  maximum  cross  range  is  achieved  with  a  bank 

(2) 

angle  of  45  deg,  as  observed  by  Slye,  '  but  only  at  low  lift-to-drag 
ratios.  At  high  L/Ds  with  a  low  terminal  velocity,  the  maximum  side  range 
appears  to  be  attained  when  the  bank  angle  is  som.'iwhat  less  than  45  deg. 


Also  at  high  lift-to-drag  ratios,  the  vehicle  proceeds  along  its  spiraling 


path  far  enough  to  achieve  large  heading  changes,  even  completely  re¬ 
versing  its  direction  of  flight  by  the  time  it  reaches  very  low  velocities 
(0.2  u/u^  «  0.3). 

It  was  assumed  earlier  that  the  vehicle  would  not  depart  too  far 
from  its  original  entry  plane.  But  our  results  indicate  that  this 


would  not  be  a  good  assumption  for  trajectories  that  are  allowed  to 


proceed  to  low  terminal  velocities  where  the  vehicle  has  a  high  L/D 
and  employs  a  medium  bank  angle  of  '^5  deg.  Thus,  those  results  shown 


that  involve  crossrange  distances  of  more  than-2000  n  mi  or  so  are 
probably  approximate  at  best.  Hov'ever,  if  the  primary  application  of 
these  results  is  for  studying  synergetic  plane  changing,  such  a  re¬ 
striction  would  not  place  serious  limits  on  the  usefulness  of  this 
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Downrange  distance,  x  (thousands  of  nmi) 


Fig. 7— Ground  trace  during  turn,  L/D  =2 
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Fig. 8— -Ground  trace  during  turn,  L/D  =  3 
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investi^stion ,  because  in  this  case  the  heading  change  desired  may 
occur  well  before  the  vehicle  departs  very  far  from  the  original 
orbital  and  entry  plane. 
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IV.  TRAJECTORIES  FOR  BANK  ANGIZS  OF  90  DEGREES 

In  this  section  we  will  discuss  the  special  case  of  constant  bank 
angle  turns  where  0  *  90  deg.  In  this  part  of  the  analysis,  the  path 
angle  may  not  be  small,  and  the  equations  developed  above  may  not  allow 
the  analyst  to  estimate  trajectory  characteristics  with  reasonable 
accuracy.  Here,  we  will  make  other  assumptions  that  permit  closed  fosn 
solutions  to  the  equations  of  motion.  These  solutions  should  be  useful 
in  making  estimates  of  vehicle  performance  when  the  bank  angle  is  90  deg. 
As  previously  noted,  the  assumption  that  the  path  angle  is  very  small 
during  reentry  stems  from  the  fact  that  there  is  some  lift  along  the 
local  vertical,  and  this  lift  force  tends  to  keep  the  path  angle  small 
for  a  large  portion  of  the  flight  through  the  atmosphere.  In  this  dis¬ 
cussion,  we  will  assume  that  the  path  angle  may  be  nonzero,  but  that  it 
is  constant  (d  y/dt  ■  0).  These  assumptions  are  similar  to  those  made 
in  analyzing  the  ballistic  reentry  situation. Additionally,  we 
will  assume  that  the  acceleration  due  to  gravity  along  the  flight  path 
is  small  compared  to  the  deceleration  due  to  drag,  that  the  acceleration 
due  to  gravity  is  constant  in  the  region  of  interest,  and  that  the  con¬ 
trol  variables  are  constant  throughout  the  flight  (specifically,  the 
bank  angle  is  90  deg). 

Under  these  conditions  the  equations  of  motion  become 


1  ^  .  i  ^  .  S. 

2  g  ds  g  dt  W 


1 

2 


"o  S''/” 


(12) 


il  d>v  „  ^ 
g  dt  W 


1 

2 


-3h  2  ^ 

p  e  u  C,A/W 

O  la 


(13) 
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In  finding  the  altitude-velocity  profile,  the  density  factor,  a, 
is  used  to  describe  the  distribution  of  atmospheric  density,  in  that 
o  ■  e  .  From  the  geometry 


dt 


.  da 
u  sin  Y.  ^ 


and  rearranging  to  eliminate  time  as  a  variable. 


du  ^  1 
do  2 


8  pQ  CpA 

W  sin  Y 


where  y  is  the  flight  path  angle  defined  positive  in  an  upward  direc¬ 
tion.  Integrating  from  a  very  high  altitude  where  the  atmospheric 
density  is  essentially  zero,  and  the  velocity  is  the  reentry  velocity. 


u 

r 


exp  ( 


8  P, 


2W 

3  ^  sin  ^ 


-Bh. 
e  ) 


(14) 


This  result  is  shown  in  Fig.  10  for  (W/C^A)  sin  y  ■  An  auxiliary 
■calc  indicates  the  altitude  decrement  that  sho’jld  be  used  for  other 
values  of  the  ballistic  coefficient  and  path  angle.  In  the  previous 
section,  where  there  was  lift  along  the  local  vertical,  the  altitude- 
velocity  profile  exhibited  a  singularity  at  orbit  velocity.  In  the 
case  of  a  nonzero  reentry  angle,  the  singularity  exists  at  the  re¬ 
entry  velocity  because  the  altitude  at  this  velocity  w*  ,  u  sumed  to 
be  essentially  infinite.  To  avoid  difficulties  in  using  the  infor¬ 
mation  in  a  practical  problem,  we  shall  later  use  the  term  initial 


Alti»udc  (thousands  of  ft) 
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Fig.  10— Altitude  velocity  profiles  for <^*90°  or 
purely  ballistic  reentries 
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velocity  as  that  velocity  where  the  dynamic  interaction  of  the  vehicle 
and  atmosphere  begins  to  be  of  interest,  and  where  most  of  the  assump¬ 
tions  made  in  this  investigation  become  reasonably  valid. 

The  vehicle  heading  is  found  by  combining  Equations  12  and  13  to 
eliminate  time  as  a  variable,  and  then  Integrating, 


=  (L/D)  log  u  /u 
e  r 


(15) 


This  result  is  nearly  the  same  as  that  obtained  at  lower  bank  angles 
(Eq.  4)  except  that  in  this  Instance,  the  full  force  of  lift  is  a 
applied  sideways.  Fig,  2  (p,  8)  can  be  used  to  show  the  heading  angle 
for  various  velocities  and  llft-to-drag  ratios  by  setting  sin  0  ■  1.  If 
a  large  change  In  heading  Is  desired  with  only  a  small  velocity  loss  due 
to  drag,  then  a  high  L/D  value  would  be  beneficial. 

The  path  length  is  found  by  noting  that 


-Sh 


log  u/u 

(T^V-n 

\  B  W  sin  Y  / 


and  can  be  substituted  in  equation  12  to  obtain 


du  .  ,  2  ,  /  /  n2 

“  -  -  J  sin  Y  u  log^  (u/u^) 


Integrating, 


u 


2 

logg  log^  (u^/u) 

.  _.l  , 

log  u  /u 
e  r 

®  3  sin  (-y) 

B  sin  (-y) 

log^ 

- 

. 

s  ■ 


(16) 
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Agaln»  a  tlngularlty  axlstt  at  tha  raantry  valocity,  and  to  uta  tha 
raaulta  ona  auat  cooaldar  uaa  of  an  Initial  valocity  slightly  lass 
than  raantry  valocity  to  avoid  an  undafinad  path  langth. 

Tha  tiaa  of  flight  whan  tha  bank  angla  is  90  dag  say  ba  found 
by  racalling  that 


^  fl  sin  y)  u  o,  a  ■  a*®** 

Substituting  into  this  axprassion  tha  solution  for  the  velocity-altitude 
profila  darivad  aarliar. 


u/u^  *  axp 


•°o 


-f»h 


•  a 


-bo 


Than  tha  diffarantial  aquation  to  ba  solvad  bacoaas 


.  (-Bu^  sin  v)  a 


and  intagrating,  tha  solution  is  in  tha  fora  of  an  axponantial 
intagral. 

bof 

(-B  u^  sin  y)  t  -  y  a” 

bOi 

Tha  liaits  in  this  fora  ara  soaawhat  awkward,  but  aay  ba  ravisad 
in  taras  of  initial,  raantry,  and  final  valocitias  by  noting  that 
log^  u^/u  ■  bo*  Thus, 


(-6u^siny)t 


log  u^ 


log  u 
«f 


(17) 
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(8) 

is  the  final  result  and  use  of  tables^  will  yield  numerical  values 
for  the  time  of  flight.  Several  examples  are  to  be  found  in  Fig.  11. 

The  load  factor  (L/W)  on  the  glider  during  its  turn  with  a  90 
deg  bank  angle  is  found  by  noting  that 

L/W  -  1/2  p  u^  a 

°  W 

and  substituting  for  the  dummy  variable,  a, 

L/W  -  R^B  (L/D)  sin  (-y)  (u^/u^)  log^  (u^/u)  (18) 

when  the  reentry  velocity  is  about  on  the  order  of  orbit  velocity. 

The  ratio  of  load  factor  to  the  sine  of  the  reentry  angle  is  shown 

in  Fig.  12  as  a  function  of  the  velocity.  It  can  be  shown  that  the 

maximum  value  of  the  load  factor  always  occurs  at  a  velocity  corres- 

*1/2 

ponding  to  u/u^  ■  e  ■  0.603,  just  as  in  a  purely  ballistic  reentry. 
The  trends  to  be  observed  are  that  the  normal  load  factor  increases 
with  L/D  and  reentry  angle.  Thus,  to  limit  load  factors  to  reasonable 
values  for  reentries  with  ■  90  deg,  a  very  small  reentry  angle  would  be 
desirable.  Note  that  in  the  previous  analysis  where  bank  angles  were 
less  than  90  deg,  the  load  factors  were  considerably  smaller  than 
those  shown  here.  Thus,  the  normal  load  factor  is  also  quite  sensitive 
to  the  bank  angle,  when  it  is  in  the  neighborhood  of  90  deg. 

The  crossrange  and  downrange  distances  of  travel  are  found  by  a 
method  similar  to  that  used  earlier  when  the  bank  angle  was  allowed 
to  be  less  than  90  deg.  First,  note  that 
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Fig. 11-Time  of  flight  ('#'•90') 


-30- 


dy  ds  .  dx  dt 

tin  m*  ^77  ■  ^  cot  a* 


du  du 


du  du 


and  that  by  dlffarantiatlng  tha  path  length  with  rcapact  to  velocity. 


^  _  _ 1 

du  “  u  fi  tin  V  log^  (u^/u) 

and  aubatltutlng,  the  differential  equations  to  be  solved  are 

sin  [log  (u^/u)'*^®] 
u  fl  aln  V  log^  (u^/u) 

cos  [  log^(u^/u)^^*^  j 
u  B  sin  V  log^  (u^/u) 

Using  tha  substitution, 

*  •  lo8^  (u/u^)**^*^ 

these  equations  reduce  to  the  faalllar  sine  and  cosine  Integrals 

(8) 

that  are  well  tabulated. 

S'  -(L/d^  \  .1.  V  [<  !■>*.  •  '‘<F  J  <«> 

*  -(L/D)  I  .in  V  [«<5  -  «<5 

where  u^  Is  the  final  velocity  of  Interest,  u^  Is  the  vehicle  veloc¬ 
ity  veil  above  the  atvosphere  (reentry  velocity),  and  Is  that 
Initial  velocity  (soewvhat  less  than  reentry  velocity)  where  the 
Interaction  of  the  vehicle  and  the  ataosphere  begins  to  be 


du 


du 
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apprcciable  and  the  assumptions  of  this  analysis  apply.  Thus, 

u.  <  u  . 
i  r 

Figures  13,  14,  and  15  show  the  ground  traces  during  steeply 
banked  turns  with  final  velocity  as  a  parameter.  The  initial  velocities 
are  u^  ■  0.98  u^,  0.95  u^,  and  0.9  u^  respectively.  Note  that  the 
greatest  heading  change  that  can  be  accomplished  with  the  least  amount 
of  velocity  loss  requires  the  highest  lift-to-drag  ratio  (L/D  ■  4  in 
this  case).  The  side  range  achieved  for  a  given  final  velocity  is 
almost  independent  of  L/D  for  final  velocities  greater  than  about  0.7  u^. 
Thus,  the  higher  L/Ds,  if  available,  are  used  only  to  turn  the  vehicle 
faster  when  the  bank  angle  is  90  deg.  If  large  side  ranges  are  desired, 
a  more  moderate  bank  angle  would  be  used,  and  the  results  of  the  pre¬ 
vious  section  would  be  applicable. 
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X  sin  7  (  n  mi ) 

Fig.  13— Ground  trace  for  <#>  =  90*"  (u;  =0.98u 
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( juj  u )  X,  ujs  X 


X  sin  7  (  n  mi  ) 

Fig.  14— Ground  trace  for  4>  -90°  (uj  =0. 95ur 
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X  sin  "K  (  n  mi  ) 

Fig.  15— Ground  trace  for  <#>  =  90®  ( u;  »  0. 90  u 
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V.  HEATING  TRENDS 

To  gain  an  insight  into  heating  trends,  we  will  assume  that  a 
given  percentage  of  the  total  mechanical  energy  lost  during  the  turn 
due  to  drag  is  converted  into  heat  which  is  then  transmitted  to  the 
vehicle.  Gazley^^^  has  assumed  that  the  fraction  of  total  energy 
transferred  as  heat  into  the  vehicle  can  be  represented  as  a  function 
of  altitude.  We  will  also  make  this  assumption,  but  will  not  specify 
the  exact  function,  since  many  analysts  may  have  different  functions 
that  they  consider  appropriate,  and,  in  practice,  a  detailed  study 
would  require  prior  specification  of  vehicle  geometry,  position  on  the 
vehicle,  thermal  protection  system,  etc.  Using  a  general  form  of  the 
transfer  fraction  (f) ,  one  may  write  the  instantaneous  heating  rate 
indicator  as 


dt 


f  C. 


3  -gh 
u  e 


(21) 


To  find  the  Indicated  heating  rate  in  a  sliapler  form,  we  merely  sub- 
-  Bh 

stitute  for  e  in  equation  21.  For  turns  where  (J  <  90  deg. 


^  . 

dt 


2  f  (W/A)  u  (1  - 


2  2 

) 

_ 


li  ^ 

P  COS  9 


(22) 


and  for  turns  where  0-90  deg. 


^  -  2  f  (W/A)  (sin  y)  u^  log^  (u/u^)  (23) 

In  both  cases,  the  heating  rate  at  a  given  velocity  is  directly  pro¬ 
portional  to  wing  loading.  When  the  bank  angle  is  less  than  90  deg. 
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the  heat  rate  is  larger  when  the  vertical  component  of  lift  to  drag 
is  small.  This  condition  could  arise  from  either  a  steep  bank  angle 
or  a  low  vehicle  lift-to-drag  ratio,  or  both.  When  the  bank  angle  is 

90  deg,  the  heat  rate  is  larger  when  the  reentry  path  angle  is  steeper,  ^ 

and  appears  to  be  a  much  stronger  function  of  velocity  than  in  cases  ^ 

with  moderate  bank  angles.  An  indication  of  total  heat  could  be  ob** 
tained  by  integrating  these  equations,  but  is  not  included  in  this 

study.  ^ 
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VI.  FINAL  COMMENTS 


In  this  study  we  have  presented  a  niethod  of  Ccilculating  turning 
trajectories  for  reentry  vehicles.  Because  the  method  leads  to  closed 
form  solutions,  made  possible  by  certain  assumptions,  the  analysis 
has  its  limitations.  However,  the  difference  between  a  sophisticated 
computation  and  the  results  presented  here  will  be  very  small,  so  long 
as  the  vehicle  does  not  depart  too  radically  from  the  initial  reentry 
plane.  Even  at  the  larger  side  ranges,  the  model  developed  in  this 
paper  may  be  used  to  provide  first-order  estimates  of  vehicle  per¬ 
formance. 

There  are  other  assumptions  that  would  also  permit  closed  form 
solutions.  The  assumptions  that  we  have  used  have  been  made  by  many 
other  analysts,  however,  and  the  errors  from  such  sources  are  gener¬ 
ally  quite  small.  We  assumed,  for  instance,  that  the  acceleration 
due  to  gravity  was  constant  at  its  sea  level  value  in  the  region  of 
interest,  and  also  that  the  atmospheric  density  decayed  exponentially 
with  altitude.  For  both  sections  of  this  analysis,  assumptions  were 
made  about  the  flight  path  angle,  and  the  analysis  would  only  be 
nearly  correct  when  such  assumptions  apply.  Fortunately,  if  synergetic 
plane  changing  is  of  interest,  all  of  the  assumptions  made  in  this 
study  appear  to  be  reasonably  good  ones.  A  measure  of  how  reasonable 
these  assumptions  are  can  be  determined  by  comparing  our  results  with 

various  trajectory  calculations  performed  by  numerical  methods  not 

(4) 

restricted  by  linearizing  assumptions. 

The  first  comparison  is  between  two  methods  for  determining 
heading  change  as  a  function  of  vehicle  velocity.  Fig.  16  shows  the 
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Fig.  16— Comparison  of  heading  change  calculations 

by  two  methods 
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results  of  this  analysis  for  a  90  deg  bank  angle  trajectory  (solid 
lines)  and  calculations  performed  on  a  digital  computer.  The  agreement 
between  the  two  methods  is  good,  and  thus  indicates  that  one  may  use 
hand  calculations  with  confidence  for  determining  heading  changes  for 
trajectories  if  the  bank  angle  is  90  deg  and  the  llft-to-drag  ratio 
ranges  from  1.0  to  3.32. 

The  altitude-velocity  profiles  calculated  by  our  model  and  by  a 
digital  computer  are  compared  in  Fig.  17.  The  conditions  for  the  com¬ 
parison  arc  noted  on  the  figure.  Although  there  is  general  agreement, 
the  hand  model  predicts  a  higher  altitude  as  the  vehicle  velocity  de¬ 
creases  below  about  0.8  u^.  The  reason  for  this  discrepancy,  even 
though  it  is  always  loss  than  15,000  ft  in  this  case,  is  that  the  hand 
analysis  is  based  on  the  assumption  that  the  path  angle  remains  constant. 
According  to  digital  calculations,  the  path  angle  is  becoming  larger 
(in  a  downward  direction).  Thus,  it  appears  that  at  the  lower  velocities, 
the  hand  model  might  not  be  accurate  enough  for  some  users.  Fortunately 
in  this  example,  however,  there  would  not  be  nwch  interest  in  velocities 
lower  than  those  shown,  because  the  vehicle  heading  has  undergone  a 
90  deg  change  by  the  time  its  velocity  has  decreased  to  about  half  its 
reentry  velocity.  For  higher  llft-to-drag  ratios,  even  less  velocity 
would  be  lost  before  such  large  heading  changes  would  bo  achieved,  and 
the  agreement  would  probably  be  good  throughout  the  region  of  interest. 

For  comparison  of  velocity-altitude  profiles  when  the  bank  angle  is 

(4) 

less  than  90  deg,  we  have  chosen  Shaver's  "bounded”  results'  (Fig. 

17).  Hand  calculations  are  compared  with  two  trajectories  where,  on 
the  average,  the  bank  angle  is  about  85  deg  (Fig.  18).  The  agreement 
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Fig.  17— Comparison  of  velocity -altitude  profiles,<#>  *  90° 
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Fig.  18— Comparison  of  altitude -velocity  profiles,  *  85 
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is  excellent  for  the  trajectory  where  the  heading  change  is  about  30  dog. 
For  the  case  with  a  60  deg  heading  change,  agreement  is  not  as  good 
because  the  bank  angle  varied  in  the  digital  calculations.  At  the 
higher  velocities  (about  0.9  »  the  bank  angle  varied  from  83  to  85 

deg,  and  the  digital  calculations  show  the  vehicle  beginning  to  level 
off  more  than  the  hand  model  would  indicate.  At  the  lower  velocities, 
when  the  bank  angle  was  allowed  to  approach  90  deg,  the  digital  calcu¬ 
lations  tend  to  more  closely  approximate  the  hand  model  predictions. 

The  altitude-velocity  profile  is  very  sensitive  to  bank  angle  varia¬ 
tions,  as  we  observed  earlier,  and  one  should  expect  variations  to  be 
large  if  the  bank  angle  is  modulated  in  flight.  The  largest  variation 
appears  to  be  about  20,000  ft  above  the  hand  predictions.  The  digital 
results  also  indicate  a  slight  tendency  of  the  vehicle  to  "skip,”  an 
effect  that  would  not  be  evident  in  a  first-order  hand  calculation. 

Only  second-order  analyses  such  as  those  performed  by  Loh  on  an  in- 

(9) 

plane  descent  would  give  indications  of  skipping. 

In  general,  there  appears  to  be  reasonable  agreement,  good  enough 
for  advanced  planning  purposes,  between  digital  and  hand  calculations, 
so  that  the  analyst  may  make  use  of  simple  calculations  to  estimate 
trajectory  parameters  for  glide  reentry  problems. 
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